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ANALYSIS OF AN IRREGULAR BOUNDARY LAYER 
BEHAVIOR FOR THE STEADY STATE FLOW OF A 
BOUSSINESQ FLUID 


CHRISTOS SOURDIS 


Dedicated to the memory of Professor Paul C. Fife 

Abstract. Using a perturbation approach, we make rigorous the formal bound¬ 
ary layer asymptotic analysis of Turcotte, Spence and Bau from the early eight¬ 
ies for the vertical flow of an internally heated Boussinesq fluid in a vertical 
channel with viscous dissipation and pressure work. A key point in our proof is 
to establish the non-degeneracy of a special solution of the Painleve-I transcen¬ 
dent. To this end, we relate this problem to recent studies for the ground states 
of the focusing nonlinear Schrodinger equation in an annulus. We also relate 
our result to a particular case of the well known Lazer-McKenna conjecture 
from nonlinear analysis. 


1. Introduction 


1.1. The problem. In [BUS, Turcotte, Spence and Bau considered the vertical flow 
of an internally heated Boussinesq fluid in a vertical channel with viscous dissipation 
and pressure work. Starting from the basic equations for conservation of mass, 
momentum and energy in a compressible fluid, and after making various appropriate 
assumptions, they were led to the study of the following boundary value problem 
for the steady state flow: 



( 1 . 1 ) 


where A > 0 is a parameter. More specifically, the parameter A represents the (non- 
dimensional) heat addition, u is the velocity, x is the scaled position and [—1,1] is 
the horizontal cross section of the vertical channel (we refer to m for more details). 

1.2. Formal asymptotic analysis as A —> oo. A formal asymptotic analysis 
carried out in [3jD] predicts the existence of solutions to m which, as A —> oo, 
behave roughly in the following way. They converge uniformly to \JA( 1 — x 2 ) over 
fixed compacts of (—1,1); they converge (in some sense) to 



( 1 . 2 ) 
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near x = — 1 and x = 1 respectively, where Y and Z should satisfy the following 
boundary value problem: 


2 y" = y 2 - s, s > 0, 

y(0) =0, y — —> 0 as s —> oo. 


(1.3) 


To convince the skeptical reader, let us note that, letting 

e = and v = -£=«, 

V A i/A 

problem ed is equivalent to the singular perturbation problem: 

e 2 v" = v 2 — (1 — x 2 ), x £ (—1,1), 

v (-!) = u(l) = 0, 


(1.4) 


to which one can apply standard, but non-rigorous, matching asymptotic techniques 
(see for instance [23]). In these terms, \/\ — x 1 serves as an outer solution which, 
however, has an irregular boundary layer , thus creating the need for the inner solu¬ 
tions in ED- 

The familiar reader may have already observed that, after a simple normalization, 
the differential equation in ( 11.311 is non other than the Painleve-I transcendent (see 
for example [9j Ch. 5]). Despite of this fact, the study of the limit problem ( 11 . 31 ) 
is nontrivial and, in fact, has quite a history (see Proposition 12.11 herein for more 
details). Combining the results of [TTJ TO], we know that problem ed has exactly 
two solutions: Y + which is strictly increasing; Y_ which has negative slope at the 
origin and exactly one local minimum. 

Actually, analogous formal considerations leave open the possibility of existence 
of solutions converging to —y/A( 1 — x 2 ), in some sense, as A — >■ oo, but this case 
lies beyond the scope of the present article. 


1.3. Rigorous known results. To the best of our knowledge, there are no studies 
of the problem ED that link it rigorously to the limit problem (11.31) . On the 
other hand, we were truly surprised when we realized the striking similarities that 
the former problem shares with a class of extensively studied superlinear elliptic 
problems of Ambrosetti-Prodi type and the famous Lazer-McKenna conjecture that 
accompanies them. Interestingly enough, however, both problems date to the early 
1980’s. In particular, for the simplified problem: 


—Am = |«|p — A(pi(x) in f2, 
u = 0 on dil, 


(1.5) 


where D is a smooth, bounded domain of R w , p £ ^1, if N > 3, p £ (l,oo) 

if N = 1,2, and > 0 is the principal eigenfunction of —A in D with Dirichlet 
boundary conditions, the Lazer-McKenna conjecture asserts roughly that the num¬ 
ber of solutions diverges as A —> oo (see j4[ [5] and the references therein for more 
details). Remarkably, the same was also conjectured in [30j for problem (11.11) and 
was subsequently verified by Hastings and McLeod in m via a shooting argument. 
Let us point out that solutions to ED should also develop an irregular boundary 
layer, as A —> oo, since the gradient of on is nonzero (by Hopf’s boundary 
point lemma). In [4], Dancer and Yan proved the Lazer-McKenna conjecture for 
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by constructing solutions with an arbitrary number of sharp downward spikes, 
located near the maximum points of tpi and superimposed on a positive minimizing 
solution, provided that A is sufficiently large. They also studied the asymptotic 
behavior, as A —> oo, of the mountain pass solutions to m and showed that they 
have a small steep peak near the boundary (combined with the irregular boundary 
layer). In connection with this, let us note that the aforementioned increasing so¬ 
lution Y+ of (11.31) is a minimizer of the natural associated energy, while the other 
solution Y _, which has a negative peak, is a mountain pass. Even though the irreg¬ 
ular boundary layer of the problem (LSD was treated mostly as a tangential issue 
in [I] , the author still had to study the elliptic analog of (11.31) (with exponent p in 
the nonlinearity). 

It follows readily from the analysis in [4j, which was variational in nature, that 
problem m has two even solutions u± such that 

u±(x) = V -4(1 — x 2 ) — (1 — x 2 ) -2 + A~^0 (1) , (1.6) 

uniformly over fixed compacts of (—1,1), as A —> oo (here, and throughout this 
paper, Landau’s symbol 0(1) denotes a quantity which is bounded independently 
of large A)- 

(2A)-*u± ^—1 + (2A)“5s^ —)• 1±(s) in Cioc[0, oo), as A —> oo. (1-7) 

Actually, only even solutions were considered in [30]. As may be expected, the 
minus case in the above result is considerably harder to establish and, for this 
purpose, the authors had to adapt some ideas from [B], Let us emphasize that U- 
lias to be constructed as a mountain pass in the class of even functions. Indeed, 
the associated linearization on m_, for large A, has similarities with a semiclassical 
Schrodinger operator with a double-well potential; thus, it is expected to have at 
least two unstable eigenvalues in the non-symmetric class (by the instability of it_ 
in the symmetric class and tunnelling phenomena, see [l3j and Remark 13.11 below), 
contradicting |14| should U- was a mountain pass solution in the non-symmetric 
class. The mountain pass solutions, in the general class, are expected to have the 
(re-scaled) profile of Y_ at one boundary point and that of Y + at the other, for large 
A. However, this does not seem to follow directly from the analysis in [3]. Let us 
also point out that the variational approach, used for showing the above, does not 
require any knowledge of the non-degeneracy of the solutions Y± of & that is the 
absence of bounded elements in the kernel of the associated linearizations. On the 
other hand, it is essentially the non-degeneracy of the corresponding Y + that allowed 
the authors of [3] to add sharp downward spikes on top of the corresponding minimal 
solution u_)_ by means of a finite dimensional variational reduction procedure. 

1.4. The main result. In this article, using a perturbation argument, we will give 
optimal estimates for the convergence in (IE3 and also provide the missing estimates 
in the intermediate zones that are not covered by (fL6l) and & In the process, 
we will prove the non-degeneracy of the “blow-up” profile Y1, and at the same time 
provide a new proof of the fact that CEU has only Y + and Y_ as solutions, which 
was originally shown in m- 

The following is our main result. 
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Theorem 1.1. Let Y and Z be either one of the two solutions Y + and Y1 of iw. 
There exists a solution u = uyz of im such that 

! u = (2 A)*Y ({2A)i(x + 1)) + 0(A$)(x + 1), 0 < x + 1 < (2 A)~*D, 

u = ( 2 A)%Y f(2A)*(x + 1 )) + 0(A*)(x + l)i, ( 2 A)~vD <x + l<6, 

{ u = (2A)iZ ({2A)i(1 - x)) +0(A*)(l-x), 0 < 1 - x < (2 A)~*D, 


[ u = (2A)iz((2A)v(l-xi)+0{Ai)(l 

and 

u = y/ A(1 — x 2 ) + 0(1)(1 — x 2 y 2 , xG 


x)i, (2A)~iD < 1 -x < 6, 
1 + (2 A)~*D, 1 - (2A)~*d\ , 


for some constants 0 < 5 D, unifoi'mly as A —> oo (for the above notation, see 
Subsection \l.f)\ below). 

Moreover, the following a-priori estimate holds for the associated linearized op¬ 
erator: There exist constants Ai, C > 0 such that if tp G C 2 [— 1,1] and f G C[— 1,1] 
satisfy 

{ -ip" + uyzV = f, a;e (-1,1), 

^(-i) = 0 = v>(l), 

for A G (Ai, oo), then 

IMU°°(-i,i) < ca 5 ||/||l~(-i,i)- 


Although we do not show it, the above estimates are optimal as can be easily 
verified by simple scaling arguments. 

In the case of even solutions, it turns out that uy.y, is asymptotically stable 
while uy_y _ is unstable with Morse index equal to two (see Remarks l3.lll3.3l and l4.ll 
below). In the nonsymmetric case, we can tell that uy_y + and Uy + Y- have Morse 
index one (see Remark 13.21 below). 

We expect that the above a-priori estimate for the linearized operator can allow 
to extend the usual variational reduction procedure, similarly to [3], in order to 
construct new solutions to m, having an arbitrary number of downward spikes 
near the origin (each of scale A~i and at an O ((lnA)A“3^ distance from the 

others) that are superimposed on the profile of uyz, for large A > 0. 

The proof of Theorem 1 1.1 1 carries over directly to the case of the singular pertur¬ 
bation problem 

{ £ 2 u" = F(u, x), x G (a, b), 
u(a) — 0 = u(b), 

provided that the following assumptions are met: F G C 3 (lx [a, b]) and there 
exists a u 0 G C 2 (a, b) D C[a, b\ such that u 0 (a) = 0 = u 0 (b), 

F(uq(x),x)= 0, xG[a, 6], F u (uq(x),x) > 0, x G (a,b), 


F u = 0, F x < 0, F uu > 0, F ux = 0 at (0, o), 
F u =0, F x > 0, F uu > 0, F ux = 0 at (0, b). 
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1.5. Method of proof. Our strategy is to apply a perturbation argument that 
has been used in many papers in the last years. This type of argument consists of 
three main steps: Firstly, one constructs a sufficiently good approximate solution to 
the problem, then studies the invertibility properties of the associated linearization 
about this approximation, and finally captures a true solution that is close, in 
some sense, to the approximate one by some type of fixed point argument. This 
approach, however, relies heavily on the good understanding of the corresponding 
limit problems, something which is not the case here since the non-degeneracy of 
Y_ does not seem to be known. In addition, the solutions that we expect to find are 
not localized in the conventional sense, as they should develop irregular boundary 
layers. In this regard, let us point out that an extra difficulty is that the convergence 
of Y± to the square root profile is algebraically slow (see (12.41) below). 

We are able to prove the non-degeneracy of Y_ by reducing m to the ground 
state problem for a nonlinear Schrodinger equation in the half-line with zero bound¬ 
ary conditions (see (12.21) below), and take advantage of the many studies that have 
been conducted on uniqueness and non-degeneracy issues for the latter problem (see 
on i na nr]). In fact, this also allows us to give a new proof of the uniqueness of 
Y±, as was originally conjectured in m and proven by completely different tech¬ 
niques in m- Armed with the knowledge of the non-degeneracy of the blow-up 
profiles Y± , we can deal with the difficulties related to the irregular boundary layer 
behavior by adapting the perturbative approach that was developed in the recent 
papers [20) [28] , where the corresponding blow-up problem featured the Painleve-II 
transcendent. 

Let us point out that, in contrast to the problems in the aforementioned refer¬ 
ences, the instability of Y_ (recall our discussion in Subsection 11.51) suggests that 
the solutions of ( 0 ) with this blow-up profile in one of the boundaries should also 
be unstable. Therefore, the well known method of upper and lower solutions (bar¬ 
riers), see for example [26] . should not be applicable to capture such solutions of 

O'. 

1.6. Relations with geometric singular perturbation theory. The ordinary 
differential equation in (11.41) can be written as a three-dimensional, slow-fast system 
(see tm having a one-dimensional slow manifold which undergoes saddle-node 
bifurcations at the values ±1 of the slow variable x. In light of the non-degeneracy 
of Y_ that we will prove, it seems plausible that our main result can also be proven 
by the blow-up approach to geometric singular perturbation theory (see m for a 
related problem with one turning point that involves the Painleve-II transcendent). 

1.7. Extensions. We expect that an analogous result to Theorem 11.11 holds for 
positive solutions to (11.51) (at least for p > 2), having the profile of the corresponding 
to Y + one-dimensional stable solution of the blow-up problem 

( y ss + Arn-i = \y\ p - s, (s, 9) S (0, oo) x l^ 1 , 

[ y( 0,0) = 0, y(s 1 9) — s? —> 0 as s-t oo, uniformly in 9 € R^' 1 . 
orthogonal to the boundary. 

In view of the preceding discussion and the results of the current article, the only 
other solution of the above blow-up problem for which we have some non-degeneracy 
information is T_(s) for the case p = 2. However, note that this solution has 
infinite Morse index as a solution of the above problem, and thus a corresponding 
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perturbation result should involve resonance phenomena (see [7] : 22] , and especially 
m where solutions exhibiting similar irregular layered behavior were studied). 

1.8. Outline of the paper. In Section [2] we will construct sufficiently good ap¬ 
proximate solutions to ED- This is the main section of the paper and it is where we 
will prove the non-degeneracy of Y_ (the full details will be postponed to Appendix 

In Section [3j we will study the invertibility properties of the linearization of ED 
on the constructed approximate solutions, relying heavily on the non-degeneracy of 
Y±. In Section U we will use the obtained linear estimates to perturb the approxi¬ 
mate solutions to genuine ones, and also obtain related estimates for their difference 
by various comparison arguments. Finally, in Section [5j we will combine everything 
together to prove Theorem ll.il We will close the paper with an appendix, providing 
the full details of the proof of the non-degeneracy of Y_. 

1.9. Notation. In the sequel, we will often suppress the obvious dependence on A 
of various functions and quantities. Furthermore, by c/C we will denote small/large 
generic constants, independent of A , whose value will change from line to line. The 
value of A will constantly increase so that all previous relations hold. The Landau 
symbol 0(1), A —> oo, will denote quantities that remain uniformly bounded as 
A —> oo, whereas o(l) will denote quantities that approach zero as A —> oo. 

2. Construction of an approximate solution u ap 

In this section, we will construct sufficiently good approximate solutions to the 
problem ED with the same type of behavior as the solutions that we are looking 
for. 

2.1. The inner (boundary layer) solution Ui n . In this subsection, motivated 
from the aforementioned formal analysis in ED], we will use the solutions Y± of the 
blow-up problem (11.311 to construct approximate solutions to (11.111 which, however, 
are effective only near the boundary of the interval. 

The properties of the blow-up profiles Y± that we will need for the purposes of 
this paper are contained in the following proposition. 

Proposition 2.1. The boundary value problem m has exactly two solutions Y + 
and Y-. We have that [Y+)' > 0 in [0, oo), while (T_)'(0) < 0 and Y_ has a unique 
minimum at some point in (0, oo). Moreover, the solutions Y± are non-degenerate, 
in the sense that there do not exist nontrivial bounded smooth solutions of 

ip" — Y±ip = 0 in [0, oo), (/?(0) = 0. (2-1) 

Proof. Existence of two solutions, satisfying the monotonicity properties described 
in the assertion of the proposition, has been established by Holmes and Spence 
L tS! by a shooting argument (and in |?j, via the method of upper/ lower solutions 
and variational arguments, perhaps unaware of my The authors of m also 
conjectured that these solutions were indeed the only ones. Their conjecture was 
settled, to the affirmative, by Hastings and Troy m- However, their proof was, as 
we discover ( almost 25 years later!), much more complicated than necessary, and 
relied on some four decimal point numerical calculations. A truly simple proof of the 
uniqueness result of m, which in the process implies the desired non-degeneracy 
of solutions can be given, based on a previous remark of ours from [20] (see Remark 
35 therein), as follows. It is easy to see that Y + is non-degenerate and the unique 
increasing solution of (11.31) (see for example g]). Let Y be any other solution of 
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OD, and let u = Y + — Y. By an easy calculation, and the maximum principle, we 
find that u has to be a solution of 

2u" — 2Y + (s)u + it 2 = 0, s > 0, u(s) >0, s > 0, it(0) = 0, u(s) — > 0 as s — > oo. 

( 2 . 2 ) 


It has been shown recently in [8] that the general problem 
u" + -u' — V ( s)u + u p = 0, s > a, 


(2.3) 


I u(s) >0, s > a, u(a) = 0, u(s) —> 0 as s —> oo, 

where a > 0, v > 0, p € (l,oo) and V € C 1 ([a, oo)), lias at most one solution 
provided that the auxiliary function 

U(s ) = V'(s)s 3 + 0V(s)s 2 + (/? - 2 )L, 


with 


satisfies 


a = 


2v 

p + 3’ 


/3 = (p — l)a, and L = a(v — 1 — a), 


liminf U(s) > 0, 

s —^OO 


and one of the following conditions: 

(i) : U is positive in (a, oo), 

(ii) : U(a) < 0 and U changes sign only once in (a, oo). 


Moreover, if we further assume that v > 0, it has been shown in the same reference 
that the unique solution of (12.31) is non-degenerate (if such solution exists). 

We would like to adapt the proof of the aforementioned result in order to establish 
that the solution Y+ — Y- of (12.21) is unique and non-degenerate. Comparing with 
(12.31) . in the problem at hand (12.21) we have a = 0, p = 2, v = 0, a = 0, V(s) = 
Y+(s), s > 0, and U(s) = s 3 F(.(s), s > 0. We note that, by scaling, the result 
of [Sj continues to hold when a positive constant multiplies the power nonlinearity. 
Observe that the corresponding case (i) above holds. However, our potential V = Y + 
loses its positivity at s = 0 and also becomes unbounded as s —> oo. On top of that, 
in our case v = 0 and not positive as required in [5] for showing the non-degeneracy 
of the solution. Nevertheless, as we will see, the proof of [8] can be easily adapted to 
establish uniqueness for the problem (12.21) , and with some care the same can be done 
for showing that the solution Y+ — is non-degenerate. This implies at once that 
m lias exactly the two solutions Y+ and Y_, while the desired non-degeneracy 
property of the solution Y- follows readily. In Proposition IA.ll of Appendix [A] we 
will indicate how the arguments of [8] can be adapted to provide uniqueness and 
non-degeneracy for & 

The proof of the proposition is complete. 

□ 


Remark 2.1. In |18lll9j we had previously applied the same idea, used in the proof 
of Provosition [2l\ for the study ofY_, to the problem 

y" = y 2 — s 2 , s € R; y(s) — |s| —0 as |s| —> oo, 

and showed that it has exactly two solutions, one which is stable and another one 
which is unstable. Interestingly enough, during the preparation of the current article, 
we came across the paper [15] where the same result was previously obtained by 
different techniques (which are similar to those that were subsequently used in |16j ). 
Actually, the non-degeneracy of the unstable solution to the above problem, which 
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we also proved in HI and enabled us to carry out the corresponding perturbation 
analysis, does not seem to be contained in [El- 

Remark 2.2. In m it was also shown that there are solutions to that ap¬ 
proach — y/s, instead of yfs, as s —> oo. For a thorough analysis of such solutions 
and more up to date references, we refer the interested reader to |12j . 

It is easy to show that 

Y{s) — = 0(s ~ 2 ) as s —> oo, (2.4) 

and 

Y\s) =-s~^ + 0(s~ 3 ), Y"{s) = 0(s~%) as s —> oo, (2.5) 

(see also [28) App. A]). 

Let Y, Z denote either one of Y±, we define the inner solution of (ED near 
x = — 1 as 

Uin(x) = (2A)^y ^(2A)s (x + 1)^ for 0 < s = (2A)s (x + 1) < 5(2A )*, (2.6) 

where <5 > 0 is a small constant independent of A. Similarly, close to x = 1, we 
define 

u in (x) = (2A)iz ((2A)s(l - x^J for 0 < t = {2A)i{l - a:) < 6(2A)*. (2.7) 

The effectiveness of iii n as an approximate solution can be mainly measured from 
the estimate in the following proposition. 

Proposition 2.2. We have 

2 u" n — u? n + A{ 1 — x 2 ) = 0(A)(1 — x 2 ) 2 as A —> oo, 
uniformly on [— 1, — 1 + 5] U [1 — <5, 1] . 

Proof. We will sketch the proof in the case where x £ [— 1, — 1 + 5], the other case 
can be treated identically. The desired estimate follows readily by linearizing 1 — x 2 
at x = — 1, which reads as 

1 — x 2 = 2(x + 1) — (x + l) 2 , 

and using ED- 

The proof of the proposition is complete. □ 

2.2. The modified outer solution u OU f Instead of using the outer solution 

v^out —v^o^y, 

we will use a more sophisticated approximation 

Uout = |A(1 - x 2 ) - (2 A)i [s - F 2 (s)] n s ( 1 + x) - (2 A)* [t - Z 2 (t)\ ng( 1 - a;)} 2 , 

( 2 . 8 ) 

x € [—1,1], where s,t as in (USD , ED respectively (but now defined on [0, oo)), 
and ns is a smooth cutoff function such that 

f 1 if M < d, 

n d (r) = < (2.9) 

[ 0 if |r| > 2d. 

Our motivation for the definition of u out comes from niH- However, let us 
note that formulas of related nature can be found (at the formal level) in some 
books of asymptotic analysis (see [25J Ch. 8]). 
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The main result concerning u ou t is the following. 

Proposition 2.3. We have 

2 <ut - ulnt + - X 2 ) = 0(Ai)( 1 - X 2 rK (2.10) 

uniformly on [— 1 + S~ 1 (2A)~^, 1 — 6~ 1 (2A)~'s] ! as A —»• oo. Moreover, we find that 
Uout - U in = 0{Ai)( 1 - X 2 )i , ( Uout - Uin)' = 0(Ai)(l - X 2 )i, 

(2.11) 

(Uout - Uin)" = 0(A%)( 1 - X 2 )~i, 

uniformly on [—1 + <5 _1 (2A) _ b, — 1 + <5] U [1 — <5,1 — S~ 1 (2A)~i], as A —► oo. 

Proof. If x £ [— 1 + S~ 1 (2A)~i, — 1 + <5], recalling (12.81) and (12.91) , we obtain that 


u„ 


- A( 1 - x 2 ) = (2A)v [F 2 (s) - s] = 


( 2 . 12 ) 


In the same interval, we can write 

1 


Uout(x) — (2H) 5 


-(2A)i(x + l) 2 +Y 2 ((2A)i(x + l)) 


Hence, by (12.41) and (12.61) . we get 


Uout = (2 A)ey(s) 


1- -(2A)~h 2 Y~ 2 


= (2 A)*Y{s) 1 + 0(A—s)s =u in +0(As) S 2, 


uniformly for s £ [<5 1 , <5(2A) b ], as A —> oo. Moreover, direct differentiation yields 
that 


(g Afi r 

2 


(2 A)i r 
2 


-(2 A)i (x + 1) + 2(2A) b YY'] \ (2 A )b {x + l) 2 + T 2 

-s + 2(2A)iYY'{s)] f—§(2 A)~h 2 + Y 2 {s) 


and 


(2A)l 
~ 2 


-(2A)b + 2(2H)§ (T') 2 (s) + 2(24) b Y"y(s) 

-i(2H)“B S 2 + y 2 (s) 


(2A)§ r 
4 


-s + 2(2A)iYY'(s) “ -i(2H)"B S 2 + T 2 (s) 
By (12.41) . (12.51) . and (12.61) . we have that 

ilAlLsY- 1 1 - 1(2 A)~h 2 Y~ 2 


( 2 A)iY'(a) - 


1-1(2 A)~io(s) 


= u' in + 0(As) S 2 

uniformly for s £ [<5 -1 , 6(2A) b], as A —> oo. In the same fashion, we can show that 


Uout. ~ u 'ln + O(A* )s 
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uniformly for s G [6 1 , <5(2^4) s ], as A — > oo. We point out that u" ut = (u" ut ) 1 — 
(V'out )21 w ith the obvious notation, and 

(4)r = (2 A)W" + (2 A)i(Y')*Y-i+0(Ai)s-?, 

{u^) 2 = {2A)i{Y'?Y-' + 0{A*)s-l, 

uniformly as A —> oo. 

In [— 1 + S, — 1 + 25], it follows readily from (12.411 . (12.51) . and (12.811 that 

u 2 out - ^4(1 - x 2 ) = O(Ai), v! out = O(Ai), v!' out = <D(A*)> 

uniformly as A —> oo (An easy way to see these is to note that we have u ou t > cAi 
and then differentiate twice (12.121) with righthand side multiplied by the cutoff). 
Similar estimates hold for the remaining regions of [—1,1]. 

The desired assertions of the proposition follow readily from the above relations. 

□ 


The following estimates will be useful in the sequel. 


Lemma 2.1. We have 


Uout = y/A( 1 -x 2 ) + 0 ((1 - x 2 )- 2 ) (2.13) 

uniformly on [—1 + 6 _1 (2A) _ 5, — 1 + 26] U [1 — 2d, 1 — 6 _1 (2A) _ s], and 

u out = y/A(l - x 2 ), xG [-1 + 2(5,1-25]. (2.14) 

Proof. lixG [—l+6 _1 (2A) _ s, —1 + 2(5], which implies that s = (2A)*(x+l) > 6 _1 , 
from (12.81) . via (12.41) . we obtain that 


U'OUt 


A*( 1 



l + 0(«-S) 


A5( 1 



l + 0( S -i) , 


and (12.131) follows readily. Analogously we treat the case where x G [1 — 26,1 — 
6 _1 (2A) _ 5], Relation (12.141) follows immediately from the definitions (12.81) and 

m • 

The proof of the lemma is complete. □ 


2.3. Gluing the inner and outer approximations in order to create the 
global approximation u ap • We define our global approximate solution to CD 
to be the smooth function 


€ [-1, -1 + S~ 1 {2A)-i} U [1 - S~\2A)-i, 1], 


— 


Uout + (x- + X+)(uin - Uout), ze[-l + 6 : (2A) 5,1-6 1 (2A) s], 

(2.15) 

with Ui n , Uout as in (12.61) (12.71) . (12.81) respectively, and 


Xt( x ) = n S~ 1 ((2A) 5 (1 ± x)) , 


(2.16) 


where n$-i defined through (12.91) . 

The main result concerning u ap is the following. 


Proposition 2.4. Letting 


£ = 2u'f p 


U 2 ap + A{1- X 2 ), 
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we have 



0{A){ l-x 2 )\ 
0(A*)( 1-x 2 )-?, 


uniformly as A —> oo. 


x G [-1,-1 + <5- 1 (2A)-k] u [1 - S~ 1 (2A)-i, 1], 
xe [-1 + 5- 1 (2A)-s,1-(5- 1 (2A)-^], 

(2.17) 


Proof. Outside of the interpolation region [— 1+5~ 1 (2A)~^, — 1 + 2<5 _1 (2A) _ ^]U[1 — 
2S~ 1 (2A)~'s , 1 — S~ 1 (2A)~ 5 ], relation (12.171) follows at once from the assertions of 
Propositions 12.21 and 12.31 In [— 1 + 5~ 1 (2A)~i , — 1 + 26~ 1 (2A)~i], by the estimates 
of Proposition 12.21 and 12.31 we have that 

2 u" p - u 2 ap + A(1 - x 2 ) = 2v!' out - u% ut + A{1 - x 2 ) 

T U ou t ) -p 4 \_{ui n U ou t) T 2\-(Uin U 0 ut 

- 2 U out X-{Uin - Uout ) - X-(Uin ~ U out ) 2 


= 0(Az)(l + x) 2 

+0 (A^(\+x)iAi + Ai A^ (1 + x)i + A^{1 + 
+0 A?(l + x)i + A(l + x) 3 ^j 


= 0{A*){l + x) 2 , 

uniformly as A —> oo. Analogous estimates hold true in the interpolation 
[1 - 26-\2A)-i, 1 - S~ 1 (2A)-i}. 

The proof of the proposition is complete. 

We have the following two easy corollaries. 

Corrolarry 2.1. We have 


(2.18) 

region 

□ 


||£||l°°(-i,i) — II 2 u’ap — Uu P + A( 1 — a; 2 )|| i oo(_ l!l ) < CA 5 . ( 2 -19) 

Proof. It follows directly from (12.171) . □ 

Corrolarry 2.2. We have 

u ap = y/A{\ - x 2 ) + O ((1 - x 2 )- 2 ) (2.20) 

uniformly on [—1 + <5 -1 (2A)“l5, — 1 + 25} U [1 — 28, 1 — S~ 1 (2A)~i], and 

u ap = VMl-x 2 ), x G [-1 + 25,1-2(5]. (2.21) 


Proof. In view of (12.151) . if x is not in the interpolation intervals [— l+d 1 (2A) 5,-1+ 
2S~ 1 (2A)~i] and [1 — 2<5 _1 (2A) _ 5,1 — <5 _1 (2A) _ 5], the assertions of the corollary 
follow directly directly from the corresponding ones of Lemma 12.11 For x in the 
interpolation intervals, we also have to use (12.111) . 

The proof of the corollary is complete. □ 


3. Linear analysis 

In this section, we will study the linearization of m about the approximate 
solution u ap , namely the linear Schrodinger operator 

C(ip) =-ip" + Uapip, tp G C 2 (-l,l) D C[-l,l], vK±i)=o, 

(for convenience, we have divided by two). 


(3.1) 
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3.1. Properties of the potential of the Schrodinger operator. In view of 
(EU, we have 

(2 A)-iu in (-1 + (2 A)~h^j = Y(s), 0 < s < 2(2 A)i. 

Moreover, it follows from (12.111) that 

(2 A) ® (ilout — u in) 1 + (2A) 5 sj = 0(A 5 )s 5 , 
uniformly on S~ 1 ,S(2A)i , as A —> oo. Hence, via (12.151) . we find that 

(2A)~iu ap 1 + (2A)“^s^ —> Y(s) in C) oc [0,oo) as A —> oo. (3.2) 
Analogously, we find that 


(2A)' 


^1 — (2 A) —> Z(t) in Ci oc [0, oo) as A — * oo. 


(3.3) 


The asymptotic behavior of Y± (recall (11.31) ) and the definitions (12.61) (12.71) imply 
that there exist constants c,D > 0, independent of A, 5, such that 

u in > cA%( 1 - x 2 )^ 


for x £ 


— 1 + D(2A)~5, —1 + 2S~ 1 (2A)~'s u 1 — 25~ 1 (2A)-5 ,1 — D(2A)~~s 


. Ob¬ 


serve also that if x £ 

<5~ 1 ), thanks to EH and EH , we have 

-2 AH 2 


— 1 + 5 1 (2A) 5,1 — 5 1 (2 A) s (which implies that s, t > 

A. 


u~ out > A{ 1 - x 2 ) - CAUi > A{ 1 - x 2 ) - CA6*( 1 - x 2 ) > ^(1 - x 2 ) 


for some constant C > 0 independent of both A and 5, where we have decreased the 
value of 5 if necessary. Combining the above two relations with (12.111) and (12.151) . 
we arrive at 


IA 


> c\j — \/l — x 2 


X £ 


— 1 + D(2A)~s, 1 — D(2A)“5 


3.2. Uniform a-priori estimates. Let 


IMIo = IM|l°°(- 1 , 1 ) = sup \y(x)\. 

(- 1 , 1 ) 


(3.4) 

(3.5) 


The main result of this section is the following. 

Proposition 3.1. There exist constants Aq , C > 0 such that, given f £ C[— 1,1], 
there exists a unique classical solution to the boundary value problem 

C{ip) = f in (-1,1), +(±1)=0, (3.6) 

and this solution satisfies 

IMIo <CA~i\\f \\o, (3.7) 

provided that A > Aq . 

Proof. To establish existence and uniqueness for (13.61) , it suffices to show the a-priori 
estimate (13.71) which implies that the kernel of C is empty (see for example pH]). 
Suppose that the latter estimate does not hold. Then, there would exist sequences 
A n > 0, (p n £ C 2 [— 1,1], f n £ C[— 1,1] such that 

£(<Pn) = fn in (-1,1), +n(±l) = 0, 

_ 2 

A n -+ oo, Mo = 1, and A n 5 1|/„|| 0 -+ 0. 


(3.8) 

(3.9) 
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Without loss of generality, we may assume that there are x n £ (—1,1) such that 
<Pn{ x n) = ll^nllo = 1, <P' n { x n) = 0, and ip'f{x n ) < 0, 

(otherwise we can consider — <p n ). Equation (13.81) . for x = x n , gives us that 

Oap(x n ') — fn (x n ) • 

In view of (13.41) and (13.91) . we find that 


the x n ’s cannot be in ^—1 + 2<5 1 (2A n ) 5,1 — 2<5 1 (2A n ) for large n. 

Consequently, there are infinitely many n ’s such that 


£ 1,-1 + 2(5 1 (2A n ) 5 


or x n £ 


l-2<S- 1 (2A„)-5,l) . 


We may assume, without loss of generality, that the former case occurs. Therefore, 
abusing notation, we can choose a subsequence so that 


x n £ 1, —1 + 2<5 1 (2A n ) 5 


, n > 1. 


(3.10) 


Let 


<L„(s) = <Pn{x), F n (s) = fn(x ), X =-1 + (2A n ) 5.S. 
Then, relations (13.81) and (13.91) become 

- (2 A n )~iu ap (-1 + (2 An)-*s') $>n = {2A n )~iF n 
in I n = 0,2(2 A n )i , $„ = 0 on the boundary of I n , and 


(3.11) 


||$n|U~(/ n) = 1, A n *\\F n \\ L oo {In) ->0, (3.12) 

respectively. Furthermore, recalling (13.101) . we have that 

$n(sn) = 1, where s n = (2A n )i(x n + 1) £ (0,25 -1 ]. (3.13) 

Making use of (13.21) . (13.111) . (13.121) . (13.131) . and a standard diagonal compactness 
argument, passing to a further subsequence, we find that 

in Cf oc [0, oo), s n —> s* £ [0,25 1 ], 


where 

$:-y(a)$* = 0 in(0,oo), $,(0) = 0, ||$,||l-(o,oo) < 1 , and $.(«*) = 1 . 

On the other hand, by the non-degeneracy of Y± (recall Proposition 12.II) . we arrive 
at a contradiction. We have thus established the validity of (13.71) . 

The proof of the proposition is complete. □ 


Remark 3.1. Let pf < p± < ■ ■ ■ , with pf — > oo as i — > oo, denote the eigenvalues 
of the linear operators 

M±(ip) = —ip" + Y±(s)ip 

with domain {ip £ H 2 ( 0, oo), y/stp £ L 2 (0,oo), ip( 0) = 0}, which are self-adjoint in 
L 2 ( 0, oo) and have only simple eigenvalues in their spectrum since l±(s) oo as 
s — > oo (see [1]3j for more details). It follows from Prouositions \2. 1\ and \A.l\ that 

pf <0, pf > 0, while pf > 0. 
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In the case where u ap is even, using the obvious notation, we denote the eigenvalues 
of the linear operators C± in m by Xf < Xf <■■■. Arguing as in pmim , it 
follows readily that 


Af +1 -A f = 0(A~ k ) and 




O 




= 1,3,5,--- ,2 


12 


with k,n € N fixed, as A -A oo. The main observation is that, because of the 
simplicity of the eigenvalues, the associated (normalized) eigenfunction to X 2 m -i 
is even whereas that associated to X 2 m is odd, for m > 1. Thus, the eigenvalue 
problem for C in (—1,1) reduces to two eigenvalue problems in (—1,0) with boundary 
condition ip(— 1) = 0, <p(0) = 0 and y>(—1) = 0, <//(0) = 0 respectively. The main 
point being that the reduced eigenvalue problems have only one turning point (at 
x = — 1 ) and the proof of [TBJ Prop. 3.25] applies directly. We expect that, as 
in UM, ^e difference between two clustering eigenvalues is actually exponentially 
small. 


Remark 3.2. In the case where u ap is nonsymmetric, one can adapt the proof of 
Prop. 3.25] to show that the corresponding linear operator C in \3.1\ ) has only 
one negative eigenvalue, which satisfies Ai = pfA~i +0 ^ A~as A —> oo (where 

pf as in Remark \3.1\) . However, it is not clear to us how to obtain asymptotic 
expansions for the rest of the eigenvalues. Certainly this has to depend on the 
ordering between {p~} and {pf}. 

Remark 3.3. It is easy to see that relation \3. 7| j as well as the assertions of 
Remarks 1 , 9.11 and 13.21 continue to hold if the potential of C was u ap + 4> with 
A-vWWo -A 0 as A —> oo (clearly (f> has to be even for the latter remark to hold). 


4. Existence of solutions and estimates 
We seek solutions of (ED in the form 

u = u a p + 4>. (4.1) 

Substituting this ansatz in EI, and rearranging terms, we see that </> solves 
-2(t>" + 2u ap <j>=-<j> 2 + 2u'f p -u 2 ap + A(l-x 2 ), *€(-1,1); 0(± 1) = 0. (4.2) 
The next proposition is the main result of this section. 


Proposition 4.1. If A is sufficiently large, there exists a constant C > 0 and a 
unique solution of such that 


Mo < CAi. 

(4.3) 

Proof. Let us write (14.21) in the abstract form 


2C(<t>) = N{<t>) + E\ 0(±1)=O, 

(4.4) 

where C was studied in Section [3] 


M{f) = -(f 2 , 


and £ was defined in (|2.17p. 



For AI > 0, consider the closed ball of C[— 1,1] that is defined by 
B M = {f&C[- 1,1] : Mo <MAi}. 
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We will show that, if M is chosen sufficiently large, the mapping T : Bm —>■ 
C 2 [— 1,1], defined by 

C(rm= W) + £; rm± i) = o, 

(recall Proposition 13.Ill , maps Bm into itself and is a contraction with respect to 
the || • ||o norm, provided that A is sufficiently large. Let </> G Bm , via (12.191) and 
s we have 

IITWIlo < CA-§(||Af(^)||o + ||£|| 0 ) 

< CA~iM 2 Ai +CA~iAi 


< CAi(M 2 A~i + 1), 

where C > 0 is independent of both large A and M. By virtue of the above relation, 
we can choose a large M > 0 such that T maps Bm into itself, for all sufficiently 
large A. From now on, we fix such an M. Similarly, for <j>i, <f >2 E Bm , we have 

\\T(</>i) - T(</> 2 )||o < CA *is ||— 0oIIo = CA =||0i — </>o||o, 

which implies that, for large A 1 the mapping T : Bm —> Bm is a contraction. Hence, 
by Banach’s fixed point theorem, we infer that T has a unique fixed point in the 
closed set Bm- In turn, this furnishes a solution of (14.21) which satisfies the uniform 
estimate (TPll . 

The proof of the proposition is complete. □ 


Remark 4.1. Ifu ap is even, we can of course restrict ourselves to even fluctuations 
(j) in 

In the next two lemmas we will show that estimate (14.31) can be improved away 
from the boundary points. 


Lemma 4.1. Let (j) be as in Proposition EM} Given L > 1, there exists a constant 
Cl > 0 such that 


mx)\<c L Ai, 


X € 


-1, -1 + (2 A)-*l] (J [l - (2 A)-iL, 1 


for all A sufficiently large. 


(4.5) 


Proof. Let 

T(s) = (2 A)~if (-1 + (2A)"*a) , 
From (14.21) . we find that 


s G 


0,2(2 A)i 


(4.6) 


-*"+(2 A)-i Uap (-1 + (2 A)~vs) T = -^2A)~H 2 +^(2A)-i£ (-1 + (2A)~is) , 

(4.7) 

for s G ^0,2(2A)s^, and ’F(O) = 0. Furthermore, from (14. 3D . and (14.61) . it follows 
that 


Ml 


L°°(o,2(2A)5) 


< c. 


(4.8) 


In turn, relations (12.191) . (13.21) . (14.71) and (14.81) imply that, given L > 1, there exists 
a constant Cl > 0 such that 


M"(s)|<C l on [0, L] 
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provided that A is sufficiently large. Consequently, it follows from (14.81) . the above 
relation, and the elementary interpolation inequality 

II^'IU”(o,l) < 2 ||ih|| i cx,( 0ii ) + ||4 / "||l=o( 0jL ) 

(keep in mind that L > 1), that 

|*'(«)| < Cl on [0,1/], 

provided that A is sufficiently large (for some possibly larger constant Cl). Now, 
the validity of estimate (14.51) on the interval — 1,—1 + (2A)~iL follows directly 

via m • Analogously we can show its validity on the interval 1 — (2 A) sL, 1 . 
The proof of the lemma is complete. □ 

Lemma 4.2. Let (f> be as in Proposition ED There exist positive constants C, D 
such that 

\<f>(x)\ < C(1 - a; 2 ) -1 , xG \-l + (2A)~^D,l-(2A)~iD , (4.9) 

provided that A is sufficiently large. 

Proof. Let 

f(x) = -A'(1 -x 2 )-\ x€ -l + (2A)-iD,l-(2A)-iD , 

with constant K > 0 to be determined, and D > 0 to be chosen larger than that in 
(i3~H) . such that if is a lower solution to (14.21) on the above interval. Differentiating 
twice gives us 

f = - K(6x 2 + 2)(1 - a; 2 ) -3 > -8^(1 - a: 2 ) -3 . 

Recalling (12.171) . (13.41) . and (14.31) . we find that 

+ 2u ap ± + £ - £ < 16K(1-x 2 )~ 3 -cAiK(l-x 2 )-^+CAi 

+CA% (1 — x 2 )-^ 


< (1 — a; 2 ) 2 CK(1 — x 2 ) 2 — cKA^+CA 

< (1 — a; 2 ) _ 5 


CKD-iA 5 - cKAi + CA 5 


< 


(1 — a; 2 ) 2 


-\cKA* +CAi 


where the constants c, C are independent of both A and D , having increased the 
value of D if necessary. Hence, we can chose a large K > 0 such that 


—2ijj" + 2Uapif + if 2 — £ < 0, x £ 


— 1 + (2A)~i D, 1 


(2 A)~iD 


provided that A is sufficiently large. By virtue of (14.21) . (14. 31) . the above equation, 
and making use of the maximum principle, we deduce that 


—C( 1 — x 2 ) 1 < <j>{x), 


X G 


-l + (2A)-iD,l 


(2 A)~*D 


for some large constant C > 0 and all large A. Analogously we can establish the 
other side of the desired estimate (14.91) . 

The proof of the lemma is complete. □ 


In summary, we have the following. 
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Proposition 4.2. There exists a solution of EJ such that 

u-u ap = 0{A^)(l-x 2 ), ie [-1,-1 + (2A)-sd 1 U [l-(2A)-5D,l 

u-u ap = 0(l)(l-x 2 )~ 1 , xG -1 + (2A)~iD,l - (2A)-$D 
for some constant D 1, uniformly as A — > oo. 


5. Proof of the main result 

From Propositions ^. 31 and 14. 21 relation (12.151) . Corollary 12.21 and Remark 13.31 we 
can infer the validity of Theorem 11.11 


Appendix A. Uniqueness and non-degeneracy of solutions for 

problem (12.21) 

Proposition A.l. Problem it 2. ‘A) has Y + — Y1 as its unique solution. Moreover, 
this solution is non-degenerate, namely there are no nontrivial bounded solutions to 

m- 


Proof. To show uniqueness, we argue by contradiction and assume that there exist 
two distinct solutions u\ and U 2 of (12.21) . As in 13 El HU, the solutions u\ and U 2 
can be chosen such that m'^O) < u' 2 { 0) and such that they intersect at most once 
in (0, oo) (this is achieved by a shooting argument, making use only of the smooth 
dependence on the initial data of solutions to the ordinary differential equation in 
(E3). Under this assumption, as in [313 HZ], we have that 


d_ ( ui(s) \ 

ds \ U2 (s) J 


> 0 , 


s > 0. 


(A.l) 


We point out that, in the above calculation, the terms involving Y + cancel each 
other, and thus the form of Y+ is irrelevant for this part of the proof. Furthermore, 
if we define 

E(s; u ) = [u'(s)] 2 — F+(s)u 2 (s) + iu 3 (s), s > 0, u G C 2 ([0, oo)), 

o 

a direct calculation yields that 

4-U(s; Ui ) = -Ylu 2 <0, s > 0, t = 1,2. (A.2) 

ds 

Therefore, using the standard fact that any solution of (12.21) decays super-exponentially 
as s —> oo, we obtain that 


E(s\Ui) > lim E(s;Ui ) = 0, s > 0, i = 1,2. (A.3) 

s—too 

Next, as in m , we set 

F(s) = E(s; u 2 ) — ^ E{s-,u i), s > 0. 

Note that, thanks to l’hospital’s rule, we have F( 0) = 0. A direct calculation, 
making use of (IA.2I) . yields that 



E(s\ Mi), s > 0. 
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So, in view of (IA.1I) and (IA.3I) . we get that 


F'(s) >0, s > 0. 

Consequently, noting that (1A.1I) implies that 


Ms) u' 2 {0) 

«i (s) ui(0)’ 


s > 0, 


and making once more use of the super-exponential decay of U\ and u 2 , we arrive 
at the strict inequality 

0 = F(0) < lim F(s) = 0, 

s —>oo 

which is a contradiction. Hence, problem (12.21) has Y + — Y_ as its only solution. 

With some care, the non-degeneracy of ip can also be shown as in [B] (see also 
[3], [17]). The fact that V(s) = Y+(s) —> oo, as s —> oo, poses an obstruction in 
adapting some proofs of [5] to our setting (especially the second part of the proof 
of Proposition 3.1 therein). Nevertheless, the fact that Y\_ is positive on [0,oo) 
and decays to zero at an algebraic rate, see m , will allow us to bypass some of 
the arguments in [S], and in fact provide a more direct proof as follows. Firstly, 
motivated from [3], we define 

M = (r°w + Y + m ds ) 2 ’ 


and let X be the completion of Cq°(0, oo) with respect to || • ||. We note that 

/•OO 

H \\ 2 > Mi / tfds, 

Jo 

where Mi > 0 is the principal eigenvalue of 

—ip” + Y + (s)ip = pip, s > 0, ^>(0) = 0 , ip G T 2 (K). 


Let ip = Y + — Y_ be the unique solution of (12.21) . then <p is a critical point of the 
functional 

I ( u ) = + Y+{s)u 2 - ds, 

where / : X —> R is of class C 2 (here u+ = max{it, 0}). This functional has the 
mountain pass structure (see for instance [T]), and the unique solution ip of (12.211 
corresponds to a mountain pass solution. We point out that, even though the 
interval (0,oo) is infinite, compactness is restored by the property that Y+(s) —> oo 
as s —> oo (see [25]). We define the Morse index of <p as 

i(I , (p) = max{dimiL : H C X is a subspace such that 
I"(<p)(h, h) < 0 for all h € H\{ 0}} . 


It follows from the general theorem in na that 

i(I,<p)< 1- (A.4) 

In fact, for the specific equation, this can be shown in an elementary way (see [2]). 
As in MHZ], we introduce a perturbed functional 

J s (u) = I(u) - 6 J - i<M(s)u 2 ^ ds, uG X, 
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for small S > 0. By the maximum principle, we see that non-trivial critical points 
of Js are solutions to the problem 

( 2 u" — (2Y+(s) + S<fi(s)) u + (1 + 6)u 2 = 0, s > 0, u(s) >0, s > 0, 

{ ' ' (A.5) 

[ u(0) = 0, lim^^oo u(s) = 0. 

Observe that p is a solution of (IA.5I) for all 6 > 0. As in 8], our primary objective 
is to apply the arguments that were used for showing uniqueness for (12.21) in order 
to infer that ip is the only solution of (1A.5I) if <5 > 0 is sufficiently small. These 
arguments can be applied almost word for word to (1A.5[) . once we show that the 
corresponding relation to (IA.2I) holds. In other words, we have to show that 

2Yf.(s) + Sp'(s) >0, s > 0, (A.6) 

for sufficiently small S > 0 (under the assumptions of [8], recall our discussion 
following (12.31) . this was not possible and the authors had to argue indirectly). To 
this end, note that we have the following rough estimates: 

Y|(s) > min jc,-s~^j, \tp'(s)\ < Ce~ a , s > 0, 

for some positive constants c, C (the former estimate holds via (12.51) . while the latter 
from the super-exponential decay of p and (12.2(1 V We therefore deduce that (I A. (ill 
is valid if <5 E (0, to), where to > 0 is the minimum value of the function 

2C -1 min |c, | e s ; s > 0. 

Consequently, if 6 > 0 is sufficiently small, the function p = Y + — Y_ is the only 
solution to (1A.5I) . 

As in mmm, in order to show that the unique solution p of m is non¬ 
degenerate, we will argue by contradiction. So, assume that tp is degenerate. In 
view of (IA.4I) . this implies that there exists a 2-dimensional subspace H C X such 
that 

I"(p)(h, h) < 0 for all h € H. 

Since 

POO 

J' s '{u)(h, h) = I"(u)(h, h) — S / (2 u+— p(s)) h 2 ds, 

Jo 

for any u € X, h E H, we have 

pOO 

(tp)(h,h) = I"(p)(h,h) — S ph 2 ds. 

Jo 

In particular, we see that 


J' 5 '(p)(h,h) < 0 for all h € H\{ 0}, 


which implies that i(Js,p) > 2. On the other hand, since Js has the mountain 
pass structure, and (IA.5I) has tp as its only solution for small S > 0, we must have 
i(Jg,tp) < 1 for small <5 > 0. We have therefore arrived at a contradiction, thus 
completing the proof of the non-degeneracy of tp. 

The proof of the proposition is complete. □ 


Acknowledgment. The research leading to these results has received fund¬ 
ing from the European Union’s Seventh Framework Programme (FP7-REGPOT — 
2009 — 1) under grant agreement n° 245749 and by the ARISTEIA (Excellence) 


















IRREGULAR BOUNDARY LAYER BEHAVIOR IN A BOUSSINESQ FLUID 


20 


programme “Analysis of discrete, kinetic and continuum models for elastic and 

viscoelastic response” of the Greek Secretariat of Research. 

References 

[1] A. Ambrosetti, and A. Malchiodi, Nonlinear analysis and semilinear elliptic problems, 
Cambridge studies in advanced mathematics 104, Cambridge university press, 2007. 

[2] H. Berestycki, and J. Wei, On least energy solutions to a semilinear elliptic equation in a 
strip , Discrete Contin. Dyn. Syst. 28 (2010), 1083-1099. 

[3] J. Byeon, and Y. Oshita, Uniqueness of standing waves for nonlinear Schrodinger equations, 
Proceedings of the Royal Society of Edinburgh 138A (2008), 975-987. 

[4] E. N. Dancer, and S. Yan, On the superlinear Lazer-McKenna conjecture , J. Differential 
Equations 210 (2005), 317-351. 

[5] E. N. Dancer, and S. Yan, On the Lazer-McKenna conjecture involving critical and super¬ 
critical exponents, Methods Appl. Anal. 15 (2008), 97-119. 

[6] M. del Pino, and P. L. Felmer, Spike-layered solutions of singularly perturbed elliptic prob¬ 
lems in a degenerate setting , Indiana Univ. Math. J. 48 (1999), 883-898. 

[7] M. del Pino, M. Kowalczyk, and J. Wei, Concentration on curves for nonlinear 
Schrodinger equations , Comm. Pure Appl. Math. 60 (2007), 113-146. 

[8] P. Felmer, S. Martinez, and K. Tanaka, Uniqueness of radially symmetric positive solu¬ 
tions for —A u-\-u = u p in an annulus , J. Differential Equations 245 (2008), 1198-1209. 

[9] A. S. Fokas, A. R. Its, A. A. Kapaev, and V. Y. Novokshenov, Painleve Transcendents: 
The Riemann-Hilbert Approach, Amer. Math. Soc., Providence, RI, 2006. 

[10] C. Gallo, and D. Pelinovsky, On the Thomas-Fermi ground state in a harmonic potential, 
Asymptot. Anal. 73 (2011), 53-96. 

[11] S. P. Hastings, and W. C. Troy, On some conjectures of Turcotte, Spence, Bau, and 
Holmes, SIAM J. Math. Anal. 20 (1989), 634-642. 

[12] S. P. Hastings, and J. B. McLeod, Classical methods in ordinary differential equations with 
applications to boundary value problems, Graduate studies in mathematics 129, American 
Mathematical Society, (2010). 

[13] P. D. Hislop, and I. M. SiGAL, Introduction to spectral theory with applications to 
Schrodinger operators, Applied mathematical sciences 113, Springer-Verlag, New York, 1996. 

[14] H. Hofer, A note on the topological degree at a critical point of mountain pass type, Proc. 
Amer. Math. Soc. 90 (1984), 309-315. 

[15] P. Holmes, On a second-order boundary value problem arising in combustion theory, Quart. 
Appl. Math. 40 (1982/83), 53-62. 

[16] P. Holmes, and D. Spence, On a Painleve-type boundary-value problem, Quart. J. Mech. 
Appl. Math. 37 (1984), 525-538. 

[17] Y. Kabeya, and K. Tanaka, Uniqueness of positive radial solutions of semilinear elliptic 
equations in ~R N and Seres non-degeneracy condition, Comm. Partial Differential Equations 
24 (1999), 563-598. 

[18] G. Karali, and C. Sourdis, Radial and bifurcating non-radial solutions for a singular per¬ 
turbation problem in the case of exchange of stabilities, Ann. Inst. H. Poincare Anal. Non 
Lineaire 29 (2012), 131-170. 

[19] G. Karali, and C. Sourdis, Resonance phenomena in a singular perturbation problem in the 
case of exchange of stabilities, Comm. Partial Differential Equations 37 (2012), 1620-1667. 

[20] G. Karali, and C. Sourdis, The ground state of a Gross-Pitaevskii energy with general po¬ 
tential in the Thomas-Fermi limit, Arch. Ration. Mech. Anal. (2015), (DOI) 10.1007/s00205- 
015-0844-3. 

[21] N. Kawano, E. Yanagida, and S. Yotsutani, Structure theorems for positive radial solutions 
to A u + K(\x\)u p = 0 in M. N , Funkcial. Ekvac. 36 (1993), 557-579. 

[22] A. Malchiodi, and M. Montenegro, Boundary concentration phenomena for a singularly 
perturbed elliptic problem, Comm. Pure Appl. Math. 55 (2002), 1507-1568. 

[23] P. D. Miller, Applied asymptotic analysis, Graduate studies in mathematics 75, American 
Mathematical Society, (2006). 

[24] S. Nakamura, A Remark on eigenvalue splittings for one-dimensional double-well Hamilto¬ 
nians, Letters in Math. Phys. 11 (1986), 337-340. 

[25] P. H. Rabinowitz, On a class of nonlinear Schrodinger equations, Z. Angew. Math. Phys. 
43 (1992), 270-291. 



IRREGULAR BOUNDARY LAYER BEHAVIOR IN A BOUSSINESQ FLUID 


21 


[26] D. H. Sattinger, Topics in stability and bifurcation theory, Lecture Notes in Math. 309, 
Springer, Heidelberg, 1973. 

[27] S. Schecter, and C. Sourdis, Heteroclinic orbits in slow-fast Hamiltonian systems with slow 
manifold bifurcations , J. Dyn. Diff. Equat. 22 (2010), 629-655. 

[28] C. Sourdis, and P. C. Fife, Existence of heteroclinic orbits for a corner layer problem in 
anisotropic interfaces, Adv. Differential Equations 12 (2007), 623-668. 

[29] S. K. Tin, N. Kopell, and C. K. R. T. Jones, Invariant manifolds and singularly perturbed 
boundary value problems, SIAM J. Num. Anal. 31 (1994), 1558-1576. 

[30] D. L. Turcotte, D. A. Spence, and H. H. Bau, Multiple solutions for natural convective 
flows in an internally heated, vertical channel with viscous dissipation and pressure work, 
Int. J. Heat Mass Transfer 25 (1982), 699-706. 

[31] W. Walter, Ordinary differential equations, Graduate texts in mathematics 182, Springer- 
Verlag, New York, 1998. 

Department of Mathematics and Applied Mathematics, University of Crete, GR-714 

09 Heraklion, Crete, Greece. 

E-mail address: csourdis@tem.uoc.gr 



